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ABSTRACT 

The  computational  modeling  linkage  from  the 
atomistic  to  the  continuum  scales  of  homogeneous  and 
discrete  deformation  mechanisms  (twinning,  cleavage, 
micro-cracking,  stacking  faults,  et  cetera),  nano-  and 
microstructure  and  defects,  to  dynamic  failure  processes 
in  brittle  materials  should  provide  pathways  to  designing 
improved  mechanical  performance  through  controlled 
processing.  This  paper  will  review  results  from  an  internal 
multiscale  computational  program  in  first-principles 
design  of  armor  ceramics,  with  a  focus  on  an  optically 
transparent  armor  ceramic,  polycrystalline  (-200  jum 
grain  size)  aluminum  oxynitride  (AlON)  spinel,  that  uses 
a  multiscale  modeling  approach,  and  will  overview  the 
special  computational  challenges  that  are  required  for 
linking  the  vast  spatiotemporal  scales  from  the  quantum 
to  the  continuum.  The  GRAND  CHALLENGE  is  to  be 
able  to  design  materials  atom  by  atom  and  to  predict 
performance. 


1.  INTRODUCTION 

The  computational  bridge  between  the  atomistic  and 
continuum  length  scales  is  addressed  in  a  hierarchical 
fashion  (Fig.  1)  through  development  of:  1)  a  first 
principles  unit  cell  model  to  predict  the  anisotropic  elastic 
properties  of  AlON,  2)  a  classical  molecular  dynamics 
model  through  periodic  replication  of  the  unit  cell  model 
for  the  study  of  single  crystal  slip  and  twinning  dynamics, 
3)  a  single  crystal  anisotropic  plasticity  model  to  account 
for  the  kinematics  of  crystal  slip  and  twinning 
mechanisms  in  AlON,  4)  a  mesoscopic  polycrystalline 
computational  finite  element  model  that  incorporates 
single  crystal  deformation  kinematics,  and  explicitly 
includes  microcracks  that  are  represented  on  the  grain 
boundaries  using  cohesive  interface  laws  that  allow 
investigation  of  crack  nucleation,  growth,  and 
coalescence,  and  5)  a  continuum  computational  finite 
element  model  with  the  particular  challenge  of 
development  of  algorithms  for  transitioning  microcrack 
coalescence  behavior  at  the  mesoscale  to  the  continuum. 
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Fig.  1.  A  multiscale  constitutive  and  failure  model  for 
AlON  that  illustrates  the  spatiotemporal  dependence  of 
the  scales. 

Critical  in  the  development  of  the  hierarchical 
multiscale  model  outlined  above  is  the  experimental 
validation  of  the  individual  scales.  Quantum  and 
molecular  dynamics  predictions  of  the  three  independent 
anisotropic  elastic  constants  of  AlON,  Cn,  Ci2,  C44, 
should  be  consistent  with  each  other,  and  with 
experimentally  determined  values  on  few  tens  of  micron 
size  oriented  single  crystals  of  AlON,  using  for  example, 
Brillouin  spectroscopy  or  focused-ion-beam  (FIB)/SEM 
compression  and  tension  tests.  The  chief  challenge  for 
modeling  brittle  armor  materials  is  to  correctly  account 
for  the  growth  kinetics  of  microcracks,  inelastic 
deformation  mechanisms,  and  the  influence  of  defects  in  a 
multiscale  computational  environment.  The  propagation 
of  microcracks  at  smaller  scales  will  be  “coarse-grained,” 
to  higher  scales  where  global  fracture  failure  and 
fragmentation  is  observed.  Such  coarse  graining 
algorithms  can  be  validated  through  continuum- scale 
experiments  that  measure  dynamic  crack  propagation 
speeds,  mixed-mode  failure,  and  crack  bifurcation 
phenomena  using  the  in-house  coherent  gradient  sensing 
(CGS)  and/or  high  speed  imaging  techniques.  The 
development  of  consistent  coarse-graining  algorithms  for 
fracture  and  other  localization  phenomena,  which  are 
associated  with  failure  and  loss  of  material  stability,  is  a 
relatively  unexplored  but  critical  aspect  of  this  research. 
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2.  EXPERIMENTAL  OBSERVATIONS 
2.1  Dynamic  Fracture  and  Inelastic  Deformation 

The  complex  fracture  patterns  that  are  observed  in 
polycrystalline  AlON  from  an  Edge-on-Impact  (EOI)  by  a 
spherical  steel  impactor  travelling  at  -430  m/s  reveal  the 
multiscale  nature  of  fracture  in  AlON  (Fig.  2). 


Crossed  polarizers: 


Plane  light; 


AlON  being  penetrated  during  this  test.  The  transitional 
velocity  was  measured  as  1207  m/s.  Of  course,  a  major 
challenge  is  to  predict  these  experimental  observations 
using  the  multiscale  modeling  and  simulation  approach. 


Fig.  4.  X-ray  radiograph  of  penetration  into  AlON  above 
transitional  velocity  (Lundberg  and  Andersson,  2010). 


Fig.  2.  Shadowgraphs  corresponding  to  photographs  at 
various  times  after  impact  from  a  spherical  projectile 
(Strassburger  et  al.,  2006). 


(a)  (b) 

Fig.  3.  a)  SEM  of  microcleavage  fracture  from  an  EOI 
test,  b)  HRTEM  showing  a  high  density  of  dissociated 
<11 0>  dislocations  along  {111}  planes  have  been 
identified. 

A  scanning  electron  micrograph  (SEM)  image  of 
fracture  surfaces  of  post  mortem  samples  of  AlON  from 
an  EOI  experiment  (Fig.  3(a))  illustrate  the  complex 
nano-  and  micro-cleavage  patterns  from  these  dynamic 
tests.  Chen  and  McCauley  are  also  carrying  out 
systematic  indentation  studies  on  single  crystal  AlON 
grains  to  determine  the  deformation  mechanisms  in  an 
intrinsically  confined  test.  A  high  resolution  transmission 
electron  microscopy  (HRTEM)  image  of  a  cross  section 
through  an  AlON  indentation  is  illustrated  in  Fig.  3(b). 

In  addition  to  dynamic  compressive  strength 
measurements  carried  out  in  a  Kolsky  bar  arrangement, 
impact  transitional  velocity  measurements  are  also  being 
determined.  Fig.  4  shows  a  flash  X-ray  radiograph  of 


2.2  FIB/SEM  Experiments 

We  have  developed  a  custom  test  apparatus  for  in- 
situ  scanning  electron  microscope  (SEM)  mechanical 
testing.  Our  system  utilizes  a  5-axis  piezoelectric 
positioning  system  that  enables  precise  placement  and 
alignment  of  the  micro- specimens,  either  within  the 
tensile  grip  or  near  the  compression  platen.  Loads  are 
applied  using  a  high-resolution  linear  actuator  with  a 
reported  resolution  of  -1  nm  and  measured  with  a  strain 
gage  based  S-beam  load  cell.  The  load  cells  have  a 
capacity  of  10  or  100  g  with  a  resolution  of  -0.01  g. 
Specimen  load/stress  and  crosshead  displacements  are 
measured  using  a  customized  data  acquisition  program, 
while  the  specimen  strain  is  calculated  from  SEM 
micrographs  using  an  open  source  digital  image 
correlation  script  developed  by  Eberl  and  coworkers  for 
Matlab®.  We  use  focused  ion  beam  (FIB)  machining  to 
fabricate  micro-compression  and  tension  samples  with 
minimum  dimensions  of  1-10  pm.  Compression 
specimens  typically  have  a  uniform  cross-section  with  a 
2:1  aspect  ratio  while  the  tension  samples  are  higher 
aspect  ratio  “dog-bone”  specimens.  In  Fig.  5,  we  show  the 
tensile  response  for  a  Pd-based  metallic  glass  specimen 
loaded  in  tension.  The  inset  SEM  micrograph  shows  a 
typical  tension  specimen  with  a  minimum  gage  width  of 
-2  pm  and  overall  length  of  -35  pm  loaded  in  the  grip. 
Tension  and  compression  tests  of  AlON  are  pending; 
however,  this  metallic  glass  system  has  been  used 
extensively  in  the  development  and  validation  of  this 
custom  in- situ  SEM  test  stage.  In  future  work,  we  will 
examine  the  strength  and  fracture  behavior  of  single  and 
bi-crystal  specimens. 
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Fig.  5.  Stress-strain  response  for  a  Pd-based  metallic  glass 
specimen  loaded  in  tension.  The  inset  shows  an  SEM 
micrograph  of  a  FIB-machined  dogbone  specimen. 

2.3  Coherent  Gradient  Sensing  of  Dynamic  Fracture 

Despite  its  long  history,  capturing  relevant  data  and 
creating  models  for  the  dynamic  failure  of  materials  at  the 
macro-scale  is  a  non-trivial  task.  Material  failure  in  the 
form  of  cracks  or  other  localizations  is  a  discrete  event, 
reducing  the  effectiveness  of  many  traditional  material 
characterization  testing  techniques  that  rely  on 
homogenization.  Full-field  data,  particularly  when  over¬ 
driven  cracks  may  bifurcate  or  a  failure-mode  transition 
might  occur,  is  highly  informative  for  observing 
phenomena  as  well  as  developing  models  for  failure 
events. 

The  full-field  method  of  coherent  gradient  sensing 
(CGS)  (Tippur,  1991)  was  implemented  at  the  Army 
Research  Laboratory  to  study  the  dynamic  fracture  of 
transparent  brittle  materials.  CGS  is  a  full- field  imaging 
technique  that  produces  fringes  based  on  the  out-of-plane 
strain  gradient,  making  it  insensitive  to  noise  and  capable 
of  discerning  the  tip  of  a  crack  when  the  crack-opening- 
displacement  is  small.  Employing  an  8 -watt  532-nm 
wavelength  laser  coupled  with  a  Cordin  222C- 1 6  camera 
allows  framing  rates  as  high  as  10-million  frames  per 
second  (limited  by  laser  power),  which  is  critical  in 
capturing  cracks  traveling  several  kilometers  per  second. 
Crack  bifurcation  in  poly(methyl-  methacrylate)  (PMMA) 
at  a  nominal  crack  velocity  of  745  m/s  (Umberger,  2009) 
is  shown  in  Fig.  6.  CGS  and  the  method  of  caustics  have 
allowed  observation  of  high-speed  crack  propagation  in 
polymers,  glasses,  and  ceramics,  resulting  in  the 
computation  of  the  dynamic  fracture  toughness  as  a 
function  of  crack  speed  for  mixed-mode  cracks  common 
in  ballistic  penetration.  Initial  work  with  the  transparent 
ceramic  AlON  revealed  decided  thickness  effects; 
ongoing  work  seeks  to  determine  whether  this  effect  is 
due  to  the  material  itself  or  the  testing  methodology. 


Fig.  6.  CGS  images  of  crack  bifurcation  (initiation 
indicated  by  red  arrow)  in  PMMA. 

Simulation  of  dynamic  crack  propagation  at  the 
macro-scale  currently  relies  on  the  data  obtained  from  the 
macro-scale  laboratory  tests  discussed  above  along  with 
the  computational  technique  of  the  adaptive  insertion  of 
cohesive  zones.  The  adaptive  insertion  of  cohesive  zones, 
pioneered  by  Ortiz  and  his  co-workers  (Camacho  and 
Ortiz,  1996),  allows  the  connectivity  of  the  finite  element 
mesh  to  be  changed  based  on  a  specified  failure  criterion. 
Rate-dependent  cohesive  laws  for  the  fracture  of  PMMA 
were  developed  and  implemented  into  finite  element 
codes  featuring  adaptive  insertion  based  on  the  extensive 
experimentation  documented  in  Umberger  (2009). 
Current  simulations  can  accurately  predict  crack 
propagation  speeds  and  bifurcation;  validation 
experiments  to  test  the  methodology  with  multiple  mixed¬ 
mode  cracks  in  a  complex  geometry  have  yielded 
positive,  if  mesh  dependent,  results.  Ongoing 
experimental  efforts  with  AlON  and  the  mesoscale  work 
underway  at  ARL  will  hopefully  provide  a  similar  rate- 
dependent  cohesive  law  and  illuminate  the  crack 
propagation  mechanisms  (intra-  versus  inter- granular 
fracture,  localized  plasticity,  et  cetera)  critical  in 
predicting  dynamic  failure. 

2.4  Shock  Compression  Experiments 

Shock  compression  of  AlON  has  been  reported  by 
Cazamias  et  al.  (2001)  to  15  GPa,  by  Vaughan  et  al. 
(2001)  to  21  GPa,  by  Sekine  et  al.  (2003)  from  61  to  180 
GPa,  and  by  Thornhill  et  al.(2006)  from  5-89  GPa.  These 
studies  included  the  Hugoniot  Elastic  Limit  (HEL), 
Hugoniot  measurements  to  180  GPa  and  phase  transition, 
shear  strength,  and  spall  strength.  The  densities  of  AlON 
used  in  the  investigations  of  Cazamias  et  al.  (2001)  and 
Vaughan  et  al.  (2001)  varied  between  3.51-3.59  Mg/m3. 
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The  density  of  AlON  used  by  Sekine  et  al.  (2003)  and  by 
Thornhill  et  al.  (2006)  was  3.67  Mg/m3.  Dandekar  et  al. 
(2007)  analyzed  the  results  of  these  shock  compression 
investigations  to  determine  the  shear  strength  of  AlON. 
The  results  of  these  studies  may  be  summarized  as 
follows:  (i)  The  HEL  of  AlON  varies  between  10.5  and 

12.1  GPa.  (ii)  The  shear  strength  of  AlON  at  the  HEL 
range  from  3.3  to  3.9  GPa.  AlON  continues  to  maintain 
shear  strength  between  4.4  and  6.0  GPa  above  the  HEL 
even  when  the  compression  of  AlON  indicates  a  softening 
around  16.5  GPa.  (iii)  Sekine  et  al.  (2003)  observed  a 
softening  in  compression  around  50  GPa  and  postulated  a 
phase  transition  for  the  reported  shock  compression  of 
AlON  around  130  GPa.  (iv)  Spall  strength  of  low  density 
AlON  was  reported  to  be  0.14  and  1.7  GPa  for  impact 
stress  of  9.5  and  4.8  GPa  (Cazamias  et  al.,  2001). 
Thornhill  et  al.  (2006)  failed  to  detect  any  spall  strength 
in  AlON  when  shocked  to  5.4,  and  7.8  GPa,  respectively. 
The  compressive  durations  in  two  spall  measurements 
each  of  Cazamias  et  al.  (2001)  and  Thornhill  et  al.  (2006) 
were  0.36  and  0.18  ps,  and  0.35  and  1.16  ps,  respectively. 

The  puzzling/inconclusive  results  of  the  previous 
spall  experiments,  lead  us  to  postulate  a  hypothesis  that 
the  magnitude  of  spall  strength  of  AlON  is  dependent  on 
both  shock  induced  compressive  stress  and  its  duration.  A 
suite  of  experiments  have  been  designed  to  test  this 
hypothesis.  The  selected  compression  durations  for  these 
experiments  are  0.18,  0.54,  and  1.  ps  if  needed.  Spall 
behavior  of  AlON  will  be  determined  to  10  GPa  with  the 
two  selected  compression  durations.  Results  of  two  spall 
experiments  performed  at  3.1  and  3.2  GPa  are  shown  in 
Fig.  7.  The  compressive  profiles  show  unexpected 
fluctuations.  The  reason  for  the  observed  fluctuations  is 
not  currently  understood.  Similar  fluctuations  were  also 
present  in  previous  studies  (Thornhill  et  al.,  2006). 


Fig.  7.  Free  surface  velocity  profiles  in  AlON. 


the  smaller  duration  of  compression.  The  calculated 
values  of  spall  strength  of  AlON  when  shocked  to  3. 1-3.2 
GPa  for  the  compression  durations  of  0.18  and  0.54  ps  are 
1.65  and  1.09  GPa,  respectively  i.e.,  a  33%  decrease  in 
the  spall  strength. 

3.  COMPUTATIONAL  MODELING 

3.1  Quantum  Mechanics  (QM) 

The  atomic  structure  and  elastic  properties  of 
aluminum  oxynitride  spinel  (AlON)  at  high  pressure  (up 
to  40  GPa)  have  been  calculated  from  first  principles.  We 
have  assumed  an  “ideal”  stoichiometry  of  cubic  AlON 
with  35.7  mole  %  AIN  using  the  constant  anion  model 
(McCauley,  1978),  and  assuming  an  Al  vacancy  on  the 
octahedral  site.  The  elastic  constants  were  calculated  from 
independent  strains  that  were  applied  to  a  unit  cell, 
parameterizing  the  total  energy  as  a  function  of  the  strain 
and  from  a  stress-strain  relationship.  The  methods  gave 
quite  close  values  of  elastic  constants,  which  indicate 
convergence  of  energy  of  cut-offs  and  total  energy.  The 
purpose  of  the  calculations  is  to  determine  if  the  location 
and/or  segregation  of  N  atoms  in  the  unit  cell  affects  the 
elastic  properties  of  AlON.  The  calculations  have  been 
carried  out  for  two  random  (Fig.  8(a))  to  clustered  (Fig. 
8(b))  arrangements  of  nitrogen  atoms  in  the  unit  cell. 

At  ambient  conditions  a  clustered  distribution  of  N 
atoms  has  -1  eV  per  55  atoms  higher  total  energy  than  for 
a  random  distribution  and  slightly,  but  systematically 
lower  elastic  constants.  The  pressure  dependence  of  Cn, 
C\2  and  C44  for  random  and  cluster  distributions  of  N 
atoms  was  calculated  in  the  range  of  0-40  GPa  by 
performing  six  finite  distortions  of  the  lattice  and  deriving 
the  elastic  constants  from  the  strain-stress  relationship. 
The  calculated  values  of  dCn/dP  are  in  the  range  of  4.0- 

6.2  and  for  dC44/dP  -0.8-1. 5. 


(a)  (b) 


Fig.  8.  Cluster  (a)  and  random  (b)  distribution  of  nitrogen 
atoms  over  oxygen  sites. 


The  elastic  solution  provides  the  magnitude  of 
fluctuations  during  the  compression  of  AlON  even  though 
AlON  remains  globally  elastic.  The  effect  of  duration  of 
shock  wave  compression  is  clearly  seen  in  these  wave 
profiles.  The  magnitude  of  pull  back  velocity  is  larger  for 


3.2  Molecular  Dynamics  (MD) 

As  MD  simulations  are  to  be  used  to  explore  failure 
mechanisms  for  AlON,  it  is  crucial  that  the  interaction 


potential  used  in  the  MD  simulations  be  assessed  for 
suitability  in  describing  this  system.  It  is  equally 
important  that  the  simulation  supercell  adequately  depicts 
the  physical  system.  The  interaction  potential  model  we 
are  using  is  based  on  the  Buckingham  potential 


V  (r  )  =  A  exp 


C, 


(1) 


Polarization  of  the  oxygen  atoms  is  achieved  via  the 
shell  model  of  Dick  and  Overhauser  (1958). 


3.2.1  Random  Generation  of  Initial  Configuration 

In  the  constant  anion  model,  about  one-third  of  the 
aluminum  atoms  are  located  in  tetrahedral  coordination, 
the  remaining  aluminum  atoms  are  found  in  octahedral 
sites,  and  one  site  remains  vacant.  First  principles 
calculations  (Fang  et  al.,  2001)  on  an  idealized  unit  cell 
(in  which  the  vacancy  and  location  of  the  N  atoms  are 
assumed)  have  demonstrated  that  the  most  stable  structure 
of  AlON  has  the  Al  vacancy  at  the  octahedral  site,  with 
six  O  nearest  neighbors  ( V06 ).  The  positions  of  the  five 
nitrogen  atoms  and  the  one  vacancy  in  a  5  5 -atom  unit  cell 
are  not  actually  known  nor  are  they  unique.  Anisotropic 
NPT  simulations  (in  which  angles  between  cell  vectors 
were  constrained  to  be  orthogonal)  of  small  supercells 
composed  of  repeating  unit  cells  in  which  the  locations  of 
the  vacancy  and  N  atoms  were  the  same  produced  crystals 
that  were  not  cubic. 


3.2.2  Size  Effect  of  Cubicity 
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Fig.  9.  Effect  of  system  size  on  the  cubicity  of  the 
simulation  cell  in  which  the  locations  of  vacancies  and  N 
atoms  are  randomly  selected  within  individual  unit  cells. 
Average  cell  dimensions  are  shown  for  the  (a)  small 
system  and  (b)  large  system,  and  the  average  cell  angles 
of  the  (c)  small  system,  and  (d)  large  system,  with 
standard  error  bars.  The  reference  lines  indicate  the 
overall  average  values  for  that  system  size. 


values  shown  in  the  Table  1.  For  MD  systems  consisting 
of  1485  atoms,  both  shell  and  no-shell  configurations  also 
predict  a  range  in  elastic  constants.  Voigt  averages  for  the 
QM  model  more  closely  predict  values  obtained  from 
polycrystals,  and  independent  measurements  of  the  elastic 
constants  from  single  crystals  of  AlON  would  aid  in 
validating  the  QM  and  MD  models  presented  here. 


Ten  systems  of  1485  atoms  were  generated,  as  well 
as  ten  systems  of  11880  atoms,  in  which  the  locations  of 
the  vacancies  and  N  atoms  are  randomly  selected.  Each 
system  was  simulated  in  the  NoT  ensemble,  where  both 
the  simulation  cell  lengths  (a,  b  and  c )  and  the  angles 
between  cell  vectors  (a,  (3  and  y)  are  unconstrained.  If  an 
individual  system  is  cubic,  such  a  simulation  should  result 
in  the  averaged  simulation  cell  lengths  being  equal  (a  =  b 
=  c ),  and  the  angles  orthogonal  (a  =  (3  =  y  =  90°).  Fig.  9 
shows  the  average  cell  dimensions  ((a)  and  (b))  and  cell 
angles  ((c)  and  (d))  for  the  two  system  sizes.  While  no 
simulation  cell  ended  up  being  purely  cubic,  increasing 
the  system  size  clearly  improves  the  cubicity  of  the 
system,  as  the  data  points  get  closer  to  the  overall  average 
values,  with  decreased  error. 


Table  1.  AlON  elastic  constants. 


QM 

MD 

Experiment 

Cn 

296-306 

367-394 

- 

Cl2 

151-158 

224-228 

- 

C44 

164-183 

242-262 

- 

Ev 

315-341 

430-467 

307-320 

Gv 

127-139 

173-190 

123-128 

Kv 

199-207 

273-283 

206-214 

3.3  Single  Crystal  Dynamics 

A  continuum  mechanics  model  for  behavior  of  AlON 
single  crystals  has  been  developed.  Sections  3.3.1  and 

3.3.2  summarize  model  features  and  results.  A  more 
complete  description  is  forthcoming  (Clayton,  2010). 


The  cubic  elastic  constants,  Cn,  C\2  and  C44,  of 
AlON  have  been  predicted  using  QM  and  MD  methods, 
and  from  these,  Voigt  averages  (upper  bounds)  for 
Young’s,  Ev,  shear,  Gv,  and  bulk  modulus,  KVj  have  been 
determined  and  compared  with  experimental  values  from 
polycrystals  (Graham  et  al.,  1988),  (Table  1).  The  QM 
model  contains  either  a  random  or  cluster  distribution  of 
N  atoms  (Batyrev  et  al.,  2009)  resulting  in  the  range  of 


3.3.1  Material  Model 

The  structure  of  AlON  is  presumed  spinel,  with  cubic 
symmetry.  Deformation  gradient  F  is 

F  =dx/dX  =  FEFP  ,  (2) 

where  FE  accounts  for  thermoelastic  deformation  and 
Fp  accounts  for  slip  and  twinning.  Specifically, 

FpFp~l  =  ®nj  +  YdjfiyTsj  ®mi  >  (3) 
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with  yl  the  slip  rate  on  system  /,  fj  the  volume  fraction 
of  a  material  element  that  has  undergone  twinning  on 
system  /,  s  the  slip  or  twinning  direction,  and  m  the  unit 
normal  to  the  plane  of  slip  or  twinning.  For  spinel, 
twinning  shear  yT=2'm.  Slip  and  twin  systems  (Homstra, 
1960;  Mitchell,  1999)  are  shown  in  Fig.  10. 


W0 

Fig.  11.  Normalized  axial  stress  in  adiabatic  uniaxial 
strain;  data  points  for  plate  impact  of  AlON  poly  crystals. 


Secondary  systems 
{110}  &  {100}  planes 


|(iio)| 

<110>{Il0}  slip 
6  planes  x  1  direction 


3  planes  x  2  directions 


Fig.  10.  Slip  and  twin  systems  in  spinel  ceramic  crystals. 


Nonlinear  thermoelasticity  of  single  crystals  with 
defects  is  addressed.  Helmholtz  free  energy  density  ^is 
¥  =  'F(Ee  (4) 

with  Ee  =(FETFE  -1)/  2  elastic  strain,  #  temperature, 
pL  dislocation  density,  and  7jT  twin  boundary  density. 
The  explicit  form  of  (4)  includes  anisotropic  nonlinear 
thermoelasticity,  dislocation  line  energy,  and  stacking 
fault/twin  boundary  energy  (Clayton,  2009).  Kinetic 
equations  dictate  slip  rates  and  rates  of  twin  volumes. 
Resistance  to  slip  or  twinning  increases  with  shearing  rate 
and  dislocation  density.  Strength  varies  with  temperature 
or  pressure  in  proportion  to  the  elastic  shear  modulus. 


3.3.2  Model  Calculations 


Thermomechanical  responses  of  AlON  single  crystals 
and  polycrystals  subject  to  homogeneous  field  F  are 
studied.  Loading  pertinent  to  plate  impact  is  considered: 
adiabatic  uniaxial  strain  at  a  high  rate  (105/s).  In  Fig.  11, 
predicted  axial  stress  cr  (normalized  by  isentropic  bulk 
modulus  K0  =  216  GPa)  versus  V/V0  closely  fits 
experimental  data  (Dandekar  et  al.,  2007).  Calculations 
for  polycrystals  reported  in  Fig.  11  require  incorporation 
of  nonlinear  elasticity,  self-consistent  treatment  of  voids, 
and  pore  collapse  with  pressure.  Void  collapse  results  in 
the  kink  in  the  stress-volume  curve  for  0.96>V/V0  >0.91. 
The  model  overestimates  axial  stress  in  high  pressure 
regime  V/V0<  0.88,  suggesting  the  pressure  derivative  of 
the  bulk  modulus  under  shock  compression  is  lower  than 
the  value  used  in  this  model  (Graham  et  al.,  1988)  and/or 
a  need  for  elastic  constants  of  even  higher  order. 


Because  resistances  to  deformation  mechanisms 
shown  in  Fig.  10  are  not  known  a  priori,  several 
possibilities  are  explored:  equal  strengths  of  all  systems, 
exclusive  <1 10>{  1 11}  slip,  preferred  <1 10> { 1 10}  slip, 
and  preferred  <1 12>{  1 11}  twinning.  Each  is  sufficient 
for  matching  experimental  shear  strength  data  (Dandekar 
et  al.,  2007).  Representative  slip  and  twin  resistances  rc 
are  listed  in  Table  2,  accounting  for  rate  sensitivity  at 
loading  rate  y0  =  105/s  and  normalized  by  tangent  shear 

modulus  //,  the  latter  ^130  GPa  at  the  HEL.  Tabulated 
values  are  1.4%-3.8%  of  the  shear  modulus,  providing  for 
post-HEL  polycrystal  shear  strength  ^3-4  GPa.  When 
extrapolated  to  the  quasi-static  limit,  shear  resistances  are 
comparable  to  values  for  other  spinel  (Mitchell,  1999)  and 
to  Peierls  stresses  for  full  dislocations  on  {111}  planes. 
Such  extrapolation  may  underestimate  strength  in  the 
static  limit  because  rate  sensitivity  tends  to  decrease  with 
decreasing  rate  as  drag  mechanisms  decline. 


Table  2.  Dynamic  strengths  of  slip  and  twin  systems. 


^\Model 

System^\^ 

Normalized  shear  strength  rc  /  p 

All 

Equal 

Exclusive 
{111}  slip 

Preferred 
{110}  slip 

Preferred 

twinning 

{111}  slip 

0.026 

0.019 

0.029 

0.038 

{110}  slip 

0.026 

GO 

0.014 

0.038 

{100}  slip 

0.026 

OO 

0.029 

0.038 

{111} twin 

0.026 

OO 

0.029 

0.019 

Calculations  for  single  crystals  compressed  along 
[100]  and  [110]  demonstrate  that  maximum  shear  stresses 
can  vary  by  over  a  factor  of  two  at  compression 
corresponding  to  the  polycrystal  HEL  (V/V0  ~  0.97) 
depending  on  lattice  orientation  and  choice  of  model  from 
Table  2. 

3.4  Mesoscale  Model 

The  mesoscopic  length  scale  is  treated  using  parallel, 
three-dimensional  finite  elements  with  microcracks 
explicitly  represented  on  the  grain  boundaries  using 
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cohesive  interface  laws  allowing  investigation  of  crack 
nucleation,  growth,  and  coalescence.  The  basis  of  our 
model  is  the  Lagrangian  approach  and  finite  element 
meshes  representing  computational  microstructures  are 
constructed  using  Voronoi  tessellation  and  a  Monte 
Carlo-based  grain  growth  approach  (Rollett  and  Manohar, 
2004;  Wu  and  Sullivan,  2003;  Kuprat  (2000). 

The  volume-meshed  microstructures  consist  of 
grains,  which  are  treated  as  separate  bodies  in  a  multi¬ 
body  contact/interface  algorithm.  Currently,  only 
isotropic  linear  elasticity  is  being  used  for  the  bulk  grain 
material.  In  the  future,  anisotropic  linear  elasticity  and 
crystal  plasticity,  discussed  in  Section  3.3  will  be  used 
based  on  results  of  density  functional  theory  calculations 
of  AlON  for  elastic  constants  discussed  in  Sections  3.1 
and  3.2.  Microcracks  are  modeled  using  cohesive 
interface  laws,  where  a  relationship  exists  between 
traction  and  separation  of  two  surfaces  at  an  interface. 
Numerically,  the  imposed  stress  distribution  eliminates 
the  stress  singularity  at  the  “real”  crack  tip  by  imposing 
tractions  across  a  “virtual”  crack  tip,  acting  to  close  the 
crack.  The  imposed  distribution  relates  the  crack  opening 
displacement  and  tractions  and  is  referred  to  as  a  cohesive 
law.  In  this  work  a  linear-decreasing  cohesive  law  is  used. 

For  the  50-grain  simulations,  specimens  were  loaded 
with  a  range  of  constant  strain  rates.  During  the  initial 
stages  of  building  the  three-dimensional  computational 
framework,  the  microstructures  were  loaded  so  that  the 
top  boundary  of  the  mesh  was  moved  at  a  constant 
velocity  of  0.2  m/s,  or enom  -  0.20  ml s / 0.002  m  =  100  s”1 . 
Since  a  sudden  velocity  impulse  is  imposed  on  the 
microstructure,  this  corresponds  to  a  shock  loading.  Other 
boundary  conditions  could  be  applied  as  needed,  such  as 
confinement  loading.  Fig.  12  shows  the  engineering 
stress-strain  response  for  50  grain  loading  in  tension  and 
compression.  As  expected,  when  loaded  in  tension  the 
failure  strength  is  approximately  307  MPa  (a  model 
input).  In  contrast,  the  peak  compressive  strength  is  871 
MPa  at  £nom  =  100  s_1 .  This  is  an  interesting  result  because 
Paliwal  et  al.  (2006)  measured  720  MPa  at  quasi-static 
loading  rates  ( enom  =  103  „v~' ). 


Strain  (%) 


Fig.  12.  50-grain  simulation  in  tension  and  compression. 


Fig.  13.  Stress  contour  results  in  a  50-grain  simulation. 

not  reached.  Fig.  13(b)  shows  an  image  of  the 
microstructure  at  35.64  jus  rendered  transparent  to  show 
how  most  of  the  grain  boundaries  are  still  carrying  load 
(have  high  stress,  colored  red).  This  is  also  evident  in  Fig. 
13(c),  which  shows  a  cross  section  of  the  micro  structure 
also  at  35.64  ps. 


Contour  plots  of  stress  in  the  Y-direction  (loading 
axis)  of  the  compression  simulation  at  various  times  are 
shown  in  Fig.  13.  At  27.66  jlis,  stress  concentrations  at 
grain  boundaries  and  triple  junctions  are  pronounced  until 
32.45  ps  when  a  multiple  cracks  are  activated  and  begin 
to  open.  Once  cracks  are  activated  and  begin  to  propagate, 
it  is  not  until  35  ps,  or  2.55  ps  later  that  the  compressive 
peak  strength  is  reached.  It  is  interesting  to  note  that  at 
peak  stress  (which  would  occur  somewhere  between 
34.58  and  35.64  ps  in  the  shown  contour  plots)  the 
maximum  crack  density  (all  grain  boundaries  fractured)  is 


CONCLUSIONS 

Multiscale  modeling  of  ceramics  can  provide 
valuable  tools  to  understand  the  mechanics  and  physics  of 
the  brittle  material  inelastic  response.  Experimental 
observations  at  each  of  the  spatiotemporal  scales  shown  in 
Fig.  1  can  be  used  for  validation  of  the  computational 
models  developed  at  that  scale.  More  specifically, 
interfacial  grain  boundary  properties  and  estimated  quasi¬ 
static  strengths  of  slip  and  twin  systems  (Table  2),  can  be 
quantitatively  determined  using  FIB/SEM  tests  on  AlON 
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single  and  bi-crystals.  QM  and  MD  models  were  used  to 
predict  the  cubic  elastic  constants  of  AlON,  yet  validation 
of  these  models  awaits  experiments  on  large  single 
crystals  of  AlON.  This  work  demonstrated  3-D  mesoscale 
simulations  with  fracture  explicitly  represented  using 
cohesive  zone  models;  such  models  can  be  validated 
using  CGS  methods  for  measuring  dynamic  crack 
propagation  rates  and  bifurcation  phenomena.  Although 
the  individual  models  for  the  hierarchical  multiscale 
model  for  AlON  are  still  being  combined  into  a  coherent 
framework,  this  work  established  the  feasibility  of  the 
approach  and  the  process  for  development  of  a  virtual 
failure  model  for  brittle  materials  (Fig.  12).  Future  work 
is  aimed  at  development  of  a  fully-concurrent  multiscale 
model  for  brittle  materials  for  “top-down”  computational 
design  of  armor  ceramics. 
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